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Continuous state branching processes with immigration are studied. We 



p ^ \ are particularly concerned with the associated (non-symmetric) Dirich- 

^ ' let form. After observing that gamma distributions are only reversible 

■ distributions for this class of models, we prove that every generalized 

gamma convolution is a stationary distribution of the process with suit- 
ably chosen branching mechanism and with continuous immigration. For 
such non-reversible processes, the strong sector condition is discussed in 
, terms of a characteristic called the Thorin measure. In addition, some 

Q^ \ connections with notion from noncommutative probability theory will 

' be pointed out through calculations involving the Stieltjes transform. 

od ■ 1 Introduction 

o . 

Besides its significance in the physical context, the (time-) reversibility can be thought 
of as a mathematical condition which guarantees a certain kind of 'solvability' of the 
equilibrium state and usually makes one possible to deduce explicit consequences. 
On the other hand, it is likely that the reversibility is a restrictive condition, and 
it fails for a number of stochastic models with stationary distributions of interest. 
In this paper, our attempt will be made in quantitative discussions on the degree 
of irreversibility of such systems. Let us illustrate roughly in a general setting the 
situation we will be concerned with. Suppose that we are given a Markov process 
with a stationary distribution v and generator L, say. Then consider a bilinear form 
E defined by 

^(/, 9) = -\ Lf{x)g{xMdx), f,ge D{L), (1.1) 

where D{L) is the domain of L. The Dirichlet form is a suitable extension of £ and it 
is well-known that the symmetry of S is interpreted as the reversibility of the Markov 
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process. We say that 8 satisfies the strong sector condition if there exists a finite 
constant C such that 

S{f, 9) < C£{f, fy/'Sig, gf'\ f, g e D{L). (1.2) 

Based on a weaker version called the weak sector condition, the theory of symmetric 
Dirichlet forms has been successfully extended to non-symmetric cases in [12]. The 
strong sector condition is known also to play an essential role in the proof of the 
invariance principle for additive functionals of non-symmetric Markov processes. See 
[15] and [23]. (See also [9] for a generalization.) Intuitively, the validity of this 
condition tells us that the process is a small perturbation from a symmetric one. It 
seems typical that verification of (1.2) depends heavily on the mathematical structure 
of the process. In our subsequent discussions it will be convenient to denote by Sect(£^) 
the infimum of C's satisfying (1.2) if any, and set Sect(£) = oo otherwise. We call 
Sect(£^) the sector constant of S. Clearly Sect(£^) > 1. If £ is symmetric, we have 
Sect(£^) = 1. As expected, the converse holds true in general. (Sec Proposition 3.1 
below for the proof.) Hence the difference Sect(£) — 1 can be thought of as a 'degree' 
of asymmetry of £ and of irreversibility of the process. Our objective is to show such 
a property of Sect(£^) in an explicit way for some specific class of models. We thus 
seek for an upper bound of Sect(£^) — 1 which should be given in terms of certain 
characteristics of the models, and the bound is then required to vanish precisely in 
the reversible case. 

As a model which will be discussed in the present paper, we adopt the continuous 
state branching process with immigration (called also the CBI- process) , which is a 
Markov process on R+ := [0, oo). Fundamental results of this process, including 
limit theorems from Galton- Watson processes with immigration and the complete 
determination of the generator, are obtained by Kawazu and Watanabe [7]. Since 
then, this model has been studied extensively not only because of the rich and nice 
mathematical structure which has allowed us to obtain a number of concrete results 
of interest but also of its importance in various applications. The aforementioned 
authors showed its interesting applications in the context of stochastic analysis as 
well. In addition, since the CIR model [1] , a mathematical finance model for evolution 
of interest rate, is included as a special case (in fact, the diffusion case), the class 
of CBI-processes serves also as a useful generahzation of the CIR model in such a 
context [4]. Wc intend to reveal further aspects of the CBI-process regarding the 
non-reversible stationary distribution and the non-symmetric Dirichlet form. 

The time evolution of the CBI-process in general incorporates two kinds of dy- 
namics; the one describing the branching of particles and the other being due to 
immigration. It is of essential importance to take into consideration the effect of 
immigration. One of consequences of introducing immigration is ergodicity of the 
process; it may exhibit the strong convergence to a unique stationary distribution, 
if any, as time goes to infinity. Actually, under suitable assumptions, the positivity 
of the spectral gap of L follows from the result in [14]. (See the discussion in the 
paragraph preceding to Lemma 2.1 below for the precise statement.) Furthermore, 
our process has so nice structure as to make it possible to get information of the 
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stationary distribution through an exphcit representation of the Laplace transform. 
This formula, a key tool throughout this paper, is due to Ogura [14], who carried out 
detailed calculations of the spectral representation for the CBI-process. However, it 
is typically difficult to deduce direct expressions (e.g. the density function) of the 
stationary distribution and so one needs to exploit other structures. (Among excep- 
tions are gamma distributions, which are reversible distributions of the CIR models.) 
Another feature of the model which is crucial to us is the branching property, mean- 
ing that the law of the sum of two identical and independent processes starting from 
Xi and X2 respectively coincides with the law of a process starting from xi + X2- By 
virtue of this property, the law of the process at an arbitrarily fixed time is necessarily 
infinitely divisible and so is the stationary distribution. (See [8] for recent studies of 
stationary distributions of the CBI- processes.) In [21], the Poincare inequahty for a 
class of infinitely divisible distributions on R+ (and more general spaces) was proved 
by reducing it to an analogous estimate for the associated Levy measure. It will turn 
out that a suitably modified argument works well for the sector constant estimate. 

Non-reversible stationary distributions we will focus attention on are generalized 
gamma convolutions [2] (GGC's for short), namely weak limits of finite convolutions 
of gamma distributions. (See also Section 5, Chapter VI of [22] for general accounts 
and [6] for a recent survey and related topics.) Wc regard these distributions as 
'perturbations' from gamma distributions. A GGC without 'translation term' is 
determined uniquely by the so-called Thorin measure, which appears in the logarithm 
of the Laplace transform and prescribes the weight of convolutions. For example, 
every gamma distribution has a degenerate Thorin measure. Therefore, the actual 
problems we are going to consider in the subsequent sections are outlined as follows. 

(I) Show that there does not exist a (nondegenerate) reversible distribution of the 
CBI-processes except gamma distributions. 

(II) Given a GGC with Thorin measure m, choose a branching mechanism so that 
the CBI-process has the GGC as a unique stationary distribution. 

(III) For the bilinear form S associated with that process, give an upper bound 
C = C(rri) of Sect((£') such that C{m) = 1 if and only if m is degenerate. 

We will see that the reversibility problem (I) reduces to solving certain functional 
equations involving 'characteristics' of the mechanisms of branching and immigration. 
Our solution to (II) will turn out to rely on the theory of Bernstein functions [18]. 
We also make use of Stieltjcs transforms in order to get further information (e.g., 
the one needed to solve (HI)) on the branching mechanism chosen. In this context 
some connections with notion from non-commutative probability theory (such as the 
so-called Boolean convolution and the free Poisson distribution) will be pointed out. 

The organization of this paper is as follows. In the next section a precise de- 
scription of CBI-processes is given and then the problem (I) is solved. In Section 3, 
we present some basic results on the strong sector condition for a subclass of CBI- 
processes for which an integration by parts formula is available. In Section 4, both 
the problems (II) and (III) are solved by constructing the CBI-process associated 
with a GGC and then applying the results in Section 3. In Section 5, we give some 
examples to illustrate consequences of our results and discuss related topics. 
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2 The model and its stationary distribution 



Following [7], we begin with a precise description of our model, namely the CBI- 
process in terms of the generator. For the purpose of this paper, we shall restrict 
the discussion to a class of conservative CBI-processes. In view of Theorem 1.1', 
Theorem 1.2 in [7], and results (especially Proposition 1.1) in [14], the assumptions 
made below are not optimal but useful in order that the results of this section are 
not more complicated than are necessary in the subsequent sections. (Recently, the 
detailed analysis of stationary distributions was done in [8] for conservative CBI- 
processes.) The generator L of our process takes the following form: 

roo 

Lf(x) = ax fix) - bxf'ix) + X / [fix + y) - fix) - yf'ix)] midy) 

Jo 

roo 

+Sf'ix)+ [fix + y)-fix)]n2idy), x e R+, (2.1) 
Jo 

where a > 0,b > 0,6 > 0, and measures rii and n2 on (0, oo) are supposed to satisfy 

/■oo r p 

/ min{y^y}ni(dy)+ / yn2idy) + (1 + logy)n2((iy) < oo. (2.2) 

Jo J{0,1) J[^:°°) 

This process approximates (asymptotically critical) Galton- Watson branching pro- 
cesses with immigration in large population limit. In this context dynamical mean- 
ing of the constants and measures appearing in (2.1) may be explained as follows. 
While a is the asymptotic variance of the offspring distributions associated with the 
branching mechanisms, b comes from the first order approximation to mean 1. S is 
the rate of change in mean of immigrating population, rii and n2 describe effects of 
big changes in population size which occur in 'macroscopic time scale' and are caused 
by branch (-death) and immigration, respectively. To avoid the triviality in discussing 
the equilibrium of the model, we make the assumption implying that both branching 
and immigration mechanisms are actually present. To be precise, defining for A > 

i?(A) = -aX^ -bX- (e-^y - 1 + Ay) n^idy) 

and 

FiX)=SX+ (l-e-^>2(%), 
Jo 

we assume throughout that 

neither = nor F = 0. (2.3) 

The functions R and F are called the branching mechanism and the immigration 
mechanism, respectively, and their interplay will be crucial in the ergodic behavior 
of the CBI-process. 

As in the literature on CBI-processes, a large amount of calculations below will be 
based on the Laplace transforms, which can be expressed in terms of the associated 
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\l/-semigroup, a one-parameter family {ijj{t, ■)}t>o of non-negative functions on R+ 
determined by the equation 

^(t,A)=i?(V'(t,A)), V'(0,A)=A (2.4) 

with A > being arbitrary. Let Tt be the semigroup of the CBI-process, and for 
every A > define a function fx on R+ by fx{x) = e~^^. Then by Theorem 1.1 in [7] 

TJx(x) = exp (^-x^|J(t, A) - |J F(^Ij(s, X))ds^ , t,X>0. (2.5) 

Ogura's formula (Eq.(1.12) with o; = in [14]) for a unique stationary distribution, 
say u, of this process is 

/ fxixHdx)^exp{-^X)), A>0, (2.6) 

provided the 'Laplace exponent' $ given by 

/•^ F(u) 

is finite for all A > 0. Conversely, if the CBI-process has a stationary distribution, 
then $(A) < oo for all A > and (2.6) holds. (See Lemma 2.1 below for the proof.) 

As shown in [14], the constant b plays an important role in studying crgodic 
properties of the process. For instance, under the assumptions 6 > and that both 
R and F are analytic at A = 0, the spectral representation of Theorem 3.1 in [14] 
implies in particular that 0, b, 2b, . . . form the discrete spectrum of —L. Assuming the 
finiteness of $ only, we will see below the convergence of the transition function as 
t — )■ oo. In such a case, the stationary distribution u is necessarily infinitely divisible 
for the reason mentioned in Introduction, and therefore $ is expressed uniquely in 
the form 

$(A) = gA + / (1 - e-^y)A(dy) (2.8) 
Jo 

for some q > (the 'translation term') and measure A (called the Levy measure) 
on (0, oo) such that mm{l,y}A{dy) < oo. (See e.g. §51 of [17].) Obviously q 
is interpreted as the infimum of the support of u. The condition that $(1) < oo is 
sufficient to guarantee that $(A) < oo for every A > since by integration by parts 

R{X) = -aX^ -bX-Xj^ (l - e"^^) rH{dy) < 0, A > 0, (2.9) 

where n{{dy) = ni{[y, oo))dy. Incidentally, we remark that (0, oo) 9 A >->■ —R{X)/X 
defines a Bernstein function with characteristic triplet (6, a, ni) in the terminology 
of [18] (Chap. 3, Theorem 3.2). Let Pt{x,dy) denote the transition function of the 
CBI-process. The following lemma gives basic observations concerning ergodicity and 
can be deduced from the results announced in [16]. The proof was given in [11]. (See 
Theorem 3.20 and Corollary 3.21 there.) 



5 



Lemma 2.1 (i) Assume that $(1) < oo and let v satisfy (2.6). Then, for each 
X e R+, ■) — >■ weakly as t ^ oo. 

(a) Suppose that the CBI-process has a stationary distribution, then $(1) < oo. 
(Hi) lfb>0, then $(1) < oo. 

In our discussion, a CBI-process is said to be ergodic if it has a (unique) stationary 
distribution, or equivalently $(1) < oo. The next proposition concerns not only 
the translation term q of the stationary distribution in the ergodic case but also 
the infimum, denoted by q{t,x), of the support of Pt{x,-). Put c = ni{dy) = 
J^yuiidy). 

Proposition 2.2 (i) If a > or c = oo, then q{t,x) = for any x e R+ and t > 0. 
Under the additional condition that $(1) < oo, it holds that q = 0. 
(ii) If a — and < b + c < oo, then for any x e R+ and t > 

q{t, x) = xe-*(''+^) + (l - 6-*^''+'=)) . (2.10) 

b + c ^ ' 

Suppose, in addition, that $(1) < oo. Then q = 6/{b + c). 

Proof. As for q{t,x), the calculation is based on (2.5) combined with a general fact 
that the infimum of the support of a probability measure /j, on R_|_ is identified with 
the 'low temperature limit' — limA^.oo(c^/c^A) log/ fxdfi. Thus 



q{t, x) = lim 

A— ^oo 



x^{t,\) + 1^ F'{^{s,X))^is,\)ds 



(2.11) 



We introduce an auxiliary function Ro{)^) = —R{\)/X, which is positive and increas- 
ing for any A > 0. It follows from (2.4) that for any t > and A > 

d,u du 
t — — = / . (2.12) 

By differentiating this identity in A 

d^p _ Rm,x)) _ ^(^,A) R,m,\)) 

9A ^> - R{X) - ~~X R^ ^ ' ^ ^ ^ ^ 

since ip{t,\) < A. Also, noting that Rq{u) G [i?o('?/'(t, A)), i?o(A)] for any u e 
[■0(i, A), A], one can deduce from (2.12) 

g-tRo(A) < ^(^' ^) < g-tRo(^(*,A))_ ^2.14) 

A 

(i) To prove that q{t, x) = 0, we only need to show that limA-j-oo |^(^) A) = for each 
t > 0. Indeed, noting that F' is decreasing and that ip{s, X) is decreasing in s and 
increasing in A by (2.13), we have for any A > 1 

< /J F'ii^is, A))|^(«, X)ds < F'm, 1)) ^{s, X)ds. 



By combining (2.13) with (2.14) 



< ^(t, A) < e-*«o(V'(*.^))^?^M!lM 



< 



dX^' ' - Ro{X) - ii?o(A)' 

which tends to as A — > oo since the assumption imphes that Ro{X) — > oo. Conse- 
quently ^{t,X) converges to boundedly and by virtue of (2.11) q{t,x) = 0. 

In the case where $(1) < oo, g = hm^^oo '&'(A) = liuix-^oo F{X)/{XRo{X)). It is 
easy to see that Iim;^_^oo -^(A)/A = 6. So we conclude that q = 0. 
(ii) It is obvious that the proof of (2.10) can be reduced to showing the following two 
asymptotics; as A — >■ oo 

^(s, A) ^ e-^''+''> for each s > 
oX 

and 

F'{i{){s, A)) 5 locally boundedly in s > 0. 

Observe that Rq{X) — )> 6 + c G (0, oo) by the assumption. Hence the first inequality in 
(2.14) implies that ■ijj{s, A) — )■ oo. By (2.14) again we have ■ijj{s, A)/A — )■ exp(— (6+c)s) 
and thus (2.13) proves the first asymptotics. The second one is a consequence of the 
following estimate; for any s e [0, t] 

\F'{iP{s,X))-S\= r ye-^^^^'^^yn^idy) < T ye-'^(*'^)%2(ciy). 
Jo Jo 

Therefore (2.10) has been estabhshed. 

The last part of the assertion (ii) can be shown by F{X)/X — >■ S and -Ro(A) b + c 
together. The proof of the proposition is complete. ■ 

It is worth mentioning the diffusion case, namely the case where a,b > and 
ni = = n2. Because of (2.3) we have S > 0, and by (2.7) the corresponding CBI- 
process, known also as the CIR model, has a unique stationary distribution with 
Laplace exponent 

$(A) = -log(l + ^A) = -/ {l-e-'y)- dy. 

a \ J a Jo y 

It is a gamma distribution with parameter {6/a,b/a), which has, by definition, the 
density proportional to x^^""'^ exp{—bx/a). This stationary distribution is reversible. 
In other words, the associated bilinear form (1.1) is symmetric. 

It would be natural to ask if there is any other case which admits a nondegenerate 
reversible distribution. The following theorem, the main result of this section, gives 
a negative answer to this question. 

Theorem 2.3 If the CBI-process with generator (2.1) has a nondegenerate reversible 
distribution, then the process coincides with a CIR model. 
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Before proving this theorem we prepare a simple lemma. In what follows, the notation 
(/) or {f{x)) will stand for the integral J-^^ f{x)u{dx) with respect to the stationary 
distribution u of an ergodic CBI-process. 

Lemma 2.4 Let A, // > be arbitrary, 
(i) For any x e R+ 

-Lh{x) = {R{\)x + F{\))h{x). (2.15) 
(a) If the CBI-process has a stationary distribution, then 

{{-L)h ■ U) = i?(A)($'(A + /x) - $'(A))(/,+,). (2.16) 

Proof. (2.15) is verified by direct calculations. Using it, we have 

{{-L)fx ■ f,) = i?(A)(xe-(^+^)-) + F(A)(e-(^+'^)-). (2.17) 

Also, (2.6) and (2.7) give {xe~^^+^'>) = $'(A + /x)(e-(^+^)^) and F(A) = -i?(A)$'(A), 
respectively. (2.16) follows by plugging these equalities into (2.17). ■ 

Proof of Theorem 2. 3. By the assumption, we have the symmetry of the Dirichlet 
form. In particular, {{—L)fx ■ f^) = {[—Vjf^ ■ fx) for all A,/x > 0. By virtue of 
Lemma 2.4 (ii), this becomes 

i?(A)($'(A + ^x)- $'(A)) = i?(/.)($'(A + ^x)- ^\^x)). 

Because of (2.7), the above equahty is rewritten into 

(i?(A) - R{^))F{\ + /.) = (F(A) - F{^))R{\ + /.). (2.18) 
Differentiating in A yields 

E!{\)F{\ + ^J) + (i?(A) - R{n))F'{X + //) = F'{X)R{X + /x) + (F(A) - F{pl))R'{X + n). 
By interchanging the roles of A and /x 

R'{li)F{X + + {R{ii) - R{X))F'{X + //) = F'{ii)R{X + //) + (F(/x) - F{X))R'{X + //) . 

Summing up the above two equalities, we arrive at 

(i?'(A) + R'{ii))F{X + ii) = (F'(A) + F{ii))R{X + /x). (2.19) 

Now let A 7^ /i. Then (2.18) and (2.19) together with Cauchy's mean value theorem 
imply further that for some ^ between A and /i 

F'(A) + F'{^^) Fix) - F{^j) F'(0 " ^ ' ^ 

Here it should be noted that by (2.3) 

roo 

F\u) ^5+ ye-"%2((i|/) > (2.21) 
Jo 
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is convex and 



roo 

R\u) = -2au -b- 1/(1 - e-"y)ni{dy) < 
Jo 



is convex. Therefore, (2.20) is possible only in the case where neither F' nor —R' 
is strictly convex. Consequently, both rii and n2 must vanish. So (2.21) shows that 
5 > 0, and the positivity of b is necessary for $(A) to be finite. Moreover, a must be 
positive also since otherwise $(A) = S\/b implying that the stationary distribution 
is concentrated at S/b. The proof of Theorem 2.3 is complete. ■ 

Here are concrete examples of CBI-processes with non-reversible stationary dis- 
tributions. 

Example 2.1 (i) This example is taken from Example 4.2 of [14]. Given < o; < ^ < 
1, set a — b — 5 — 0, 

a{a + l) % , / , N dy 
Mdy) = ft; r • ^TT and n2{dy) 



Then R{\) = -A^+" and F(A) = A^. Therefore, $(A) = A^-"/(/3 - a), the Laplace 
exponent of a — Q;)-stable distribution on R+. (2.16) gives 



{{-L)h ■ /.) = 1 - 



g-*(A+M)_ 



(ii) Given < a < 1 and k > 0, define a = 0, 6 = k", 5 = 1, ^2 = and 
With these choices R{X) — — A(A -I- k)" and -F(A) = A, which together lead to 



1 — a 



r(Q;)io ' 



See e.g. [5] for information of the corresponding distribution. Note that for k = 
the stationary distribution is a (1 — Q;)-stable distribution on R+, and that as a t 1, 
$(A) tends to log(l + A/k), the Laplace exponent of a gamma distribution, provided 
that K, > 0. By (2.16) we have 

Since the class of CBI-processes studied so far seems too wide for one to obtain 
further consequences which are useful for our purpose, we will be obliged to make an 
additional restriction in the subsequent sections. In this regard, it must be remarked 
that the condition that ni = makes the correspondence between (a, 6,^2, 5) and 
(g. A) in (2.8) too simple as will be seen from the general observation below. 
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Lemma 2.5 Let a,b > and suppose that a measure n2 on (0, oo) is non-zero. For 
\>0, set 



Then for each A > 



A / (1 - e-^"')n2(dy) 



au^ + bu 



-du 



^ / c-'-yrT^idy) 



au + b 



-du. 



oo 



:i-e- 



n2{dy) 
by 



(6 = 0) 

(a = 0, 6 > 0) 



$(A) = < 



The proof requires only 'Fubini calculus' and so is left to the reader. In the light of 
this lemma, we shall proceed under the additional hypothesis that n2 = 0. 



/ (1 - e'^y) / e''"/"n^(c^^) dy (a > 0, 6 > 0) 

Jo Vio / ay 



3 Estimating the sector constant 

The main subject of this section is the estimation of the Dirichlet form. As announced 
in Introduction, we now show in a general setting that Sect(£^) > 1 holds for any non- 
symmetric Dirichlet form E. (An explicit lower bound for Sect(£^) will be discussed 
at the end of this section.) 

Proposition 3.1 Suppose that the bilinear form E in (1.1) is associated with a con- 
servative Markov process with generator L and a stationary distribution v. If E is 
non- symmetric, then Sect{E) > 1. 

Proof. We may assume that Sect(£^) < oo. Equivalently, suppose that (1.2) holds 

for some C < oo. By non-symmetry there exist f,g E D{L) such that £{f,g) > 
E{g, /). This implies E{f, /) > since otherwise (1.2) leads to the contradiction that 
^{f-iO) = ^{.9-if){= 0)- It is straightforward to see that 

^U^f + t9f A_ S{f,g)-E{gJ) 
'^^t[EifJ)Eif + tgJ + tg) V" ^ifJ) ' ^ ^ 

and hence E{f, f -\-tgY > E{f, f)E{f -\-tg, f -\-tg) for i > small enough. This shows 
that Sect(£:) > 1. ■ 

We now turn to discussing the bilinear forms E associated with ergodic CBI- 
processes. The symmetric part E{f,g) :— {E{f,g) -\- E{g, f))/2 has an expression of 
the form 

^U.g) = a{xf'{x)g'{x)) + ^{x Jni{dy){f{x + y)- f{x)){g{x + y)-g{x))) 

+^{ J n2{dy){f{x + t/) - f{x)){g{x + y) - g{x))). (3.2) 
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Here and in what follows, the domain of integration may be suppressed as long as it 
is (0, oo). (3.2) can be verified by calculating r(/, g) := {L{fg) — Lf ■ g — f ■ Lg)/2 
since £{f, g) — (r(/, g)). The main task in the remainder of this section is to give an 
upper bound of Sect(£) for a class of non- reversible CBI-proccsscs. 

In the rest of the paper, we make the restriction that ^2 = and call such pro- 
cesses continuous state branching processes with continuous immigration, henceforth 
abbreviated as CBCI-processes. Thus, F{\) = 5\ with some 5 > 0. This condition 
seems crucial in the subsequent argument, in particular, in showing an integration 
by parts formula described in Proposition 3.2 below. The notation n is used instead 
of ni and thus a measure n on (0, oo) is assumed to satisfy / mm{y'^ ^ y}n{dy) < oo 
according to (2.2). In the discussion below, we shall suppose the existence of a unique 
stationary distribution and introduce a one-parameter family of the convolution semi- 
group {vs '■— T^*^}s>o with 1/ having the Laplace exponent 

$(A) = /' = qX+ f(l- e-^y)A(dy), (3.3) 

■^0 au + b+ {1- e-''y)n{dy) 

where q > and A is a Levy measure on (0, oo). Accordingly, 5$ is the Laplace 
exponent of vs, which is a unique stationary distribution of the process with generator 

Lsfix) := ax fix) - hxf{x) +xj [f{x + y) - f{x) - yf{x)] n{dy) 
+Sf'{x). 

We call it the CBCI-process with quadruplet {a,h,n,5). This subclass of CBI- 
processes is one-dimensional version of the model discussed in [20] and [21]. We em- 
phasize that an explicit formula for the Levy density dK/dy was obtained in Lemma 
2.5 of [21] under the additional hypothesis that a, 6 > and c = / n{dy) < oo. (In 
this case g = by Proposition 2.2 (i).) An analogue of that formula is available also 
in the case where a — and < 6 -|- c < oo. Indeed, by differentiating (3.3) 

$'(A) = = /e-^^n(ci|/)) , 

b + c- je-^ynidy) ^ + c ^,b + c\b + cJ ') 

which shows that q= l/{b + c) and 

1 °° 1 

A(*) = -E(^^fi-(d,). (3.4) 

The notation will stand for the integral with respect to and the associated 
bihnear form is denoted by £^ , namely £^{f, g) — {{—Ls)f • g)s. Let To be the hnear 
hull of {/a : a > 0}. As remarked after Theorem 1.1' in [7], is a core of the 
generator of the CBI-process. 
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Proposition 3.2 Let $ be defined by the first equality in (3.3) and suppose that 
$(1) < oo. Then for each 5 > 

£\f,g) = a{xf{x)g\x))s 



+ 



{x J n{dy)f{x + y){g{x + y) - g{x)))s, f,ge J'o- (3.5) 

Proof. It suffices to sfiow (3.5) for f — f^^ g — f^j^ witfi A, > being arbitrary. We 
begin witfi a version of (2.16) in Lemma 2.4: 

{{-Ls)fx ■ f,)s = 5$'(A + /.)(W). (r{X) + + ' ^^-^^ 
wfiere R{X) = -aA^ - b\ -XJn{dy){l - e'^y). Observing from (3.3) that 
^.(/^ ^) ^a{X + fx) + b + J n{dy){l - e'^^^'^)^), 

we have 

m + ^7(X^ = a\^i + \J n{dy) (e"^^ - e-(^+^)^) . (3.7) 

Since 5^'{\ + ii){fx+^)5 = {xfx+u.{x))5, we get by plugging (3.7) into (3.6) 

{{-L5)fx ■ f^)s = {xfx+,{x))s (aX^ + A | n{dy) (e"^^ - e-(^+^)^)) (3.8) 

= a(x/A'(x)//(x))5 + n(cit/)(-A)e-^(^+^') ^^-^i^+y) _ e-'^-))^. 

This coincides with the right side of (3.5) with fx and f^ in place of / and g, respec- 
tively. ■ 

The integration by parts formula (3.5) would be interesting in its own right and 
applicable in other contexts. We here use it for the purpose of the sector constant 
estimate. While it is obvious from (3.2) that the first term on the right side of (3.5) is 
dominated by £^{f,f)^^'^£^{g,g)^^'^, the main difficulty in handling the second term 
comes from the fact that we have few information on the distribution function (or the 
density function) of the stationary distribution u^. We overcome this by a strategy 
similar to that taken in [21] for the proof of a Poincare type inequahty. That is, we 
show first that an estimate we want for the second term to satisfy can reduce to an 
analogous one for the Levy measure A, and then give a sufficient condition for the 
reduced estimate to hold. So we shall be concerned with the bilinear forms Bs{S > 0) 
and Bq on J-'q x J^q defined by 

Bsif: 9)^{xj n{dy)f{x + y){g{x + y)- g{x)))s, S > 



and 



o(/,^) = / Hdx)x I n(dy)f(x + y)(g(x + y)-g(x)), 



respectively. The first step will be done in the next theorem, the main technical 
result of this section. 
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Theorem 3.3 Let $ be defined by the first equality in (3.3). Suppose that $(1) < oo 
and that a Levy measure A on (0, oo) satisfies (3.3) with q = 0. Then, for any fixed 
< C < oo, the following two conditions are equivalent to each other: 
(i) For all5>0 

Bs{f,gf < C{xf{x)')s{x J n{dy){g{x + y)-g{x)f)s, f,g G J'o- (3.9) 
(a) For all f,g & 

Bo{f,gf <CJ A{dx)xf\xf ■ J A{dx)x J n{dy){g{x + y)-g{x)f. (3.10) 

//, in addition, a > and (3.9) holds for some 5 > 0, then '^eci{E^) < 1 + \jlC ja. 
Proof. The imphcation (i)=^(ii) foUows immediately by observing that as 5 

S~'^{xhi{x))s^ j zA{dz){hi{x + z))s^ j zA{dz)hi{z) (3.11) 

for each i e {1, 2, 3}, where hi{x) = J n{dy)f'{x + y){g{x + y) — g{x)), h2{x) = f'{x)^ 
and hslx) — J n{dy){g{x + y) — g{x)Y with f,gEJ^o being arbitrary. In (3.11) we 
have apphed the Palm formula for the underlying Poisson random measure to get 
the equality, and then used the fact that us tends weakly to the delta distribution 
at 0. (Alternatively, the equality can be verified directly for hi — fx with A > and 
extended easily. See Lemma 3.2 in [21].) 

Next, assume that (3.10) holds for all f,g ^ J^o- We must show (3.9) for every 
5 > 0. Let / and g be given as finite sums of the form f — Cifx^ and g — djfu. , 
respectively, where Cj, dj e R, Aj, > 0. In view of (3.8) 

Bsifx, f,) = {xfx+,{x))s\ J n{dy) (e"^^ - e'^'+^^y) 

= 5e-^*(^+'^)$'(A + ii)J n{dy)\ (e"^^ - e-(^+'')^) 

for any A,/x > 0. Therefore, by bilinearity Bs{f,g) equals 

5Y,Cidje-^^^^'+^^^^'(Xi + iXj) j n{dy)Xi (e-^^^ - e-(^*+'^^)^) 

= c^4e^*(^"''^) j A{dx)xe-^^*+i'i^'' j n{dy)\ (e'^'^z _ e-(A.+M,)2/j ^(3 ^2) 

where c[ = q exp(— 5<l>(Aj)), d'^ = dj exp{—d^{iJ,j)) and 

$(A, /i) = $(A) + $(/x) - $(A + /x) = j A{dz){l - e-^'){l - e"^"). A, > 0. 
Substituting the expansion 

00 rTV ^ N N 

= ^ " / A(da;i) • • • / A{dxN) - e"^'^*) - e"'^^^') 

N=0 fe=l 1=1 
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into (3.12) leads to 

oo ^Af 



Bs{f,9) = <^ E 777 / ^(^^i) • • 7 HdxN) I Hdx)x 
where 

N 

i k=l 

and 

N 

9^^,...,^Ax)^Y.d',m^-e-'^^')f,,(x) 

j /=i 

arc considered to be elements of J-q for arbitrarily given Xi, . . . ,X]sf > 0. We can 
apply now (3.10) to these functions and then use Schwarz's inequality to obtain 

oo srN 



Bsif, g) = ^Hjnl ^(^^i) • • 7 A(^^iv)So(/.„...,.,, 5.,,...,..) (3.13) 



< ^\jQi\f)\jQf\9l (3.14) 

where 

Qf\f) = <^ E ICTT / Hdxi) ■■■ A{dx^) / A{dx)xf,,_J{xy 
and Qf\g) is defined to be 

CO ^ f f f 

^^mJ ^(^^i) ■ ■ 7 HdxN) J Hdx)x j n{dy) {g 

Xl,...,XN 

{x + y)- gxi,.. .,Xn 

But analogous calculations to those for the proof of (3.13) show that 

Qf\f) = {xr{xf)s and Qf\g) ^ {x j n{dy){g{x + y) - g{x)f)s. 

Thus, (3.14) proves (3.9). 

Lastly, in view of (3.5) and (3.2) with ni — n and n2 = 0, the validity of (3.9) 
implies that 

£\f, g)<{l + j2C/^)S\f, ff's'ig, gf\ f,geJ'o 



if a > 0. This inequality is shown to extend to all functions in D{Ls), and hence 
Sect(£^) < 1 + i/2C/a as desired. The proof of Theorem 3.3 is complete. ■ 



inf{r > : n{dy) < rn{dy) in distribution sense} 



The next step is to seek conditions for (3.10) to hold. Define 

dn 
dn 

with convention that inf = oo. Clearly, this value is for n = 0. 
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Corollary 3.4 In addition to the assumptions in Theorem 3.3, assume that a > 
and that there exists a density dK/dy =: ip{y)/y such that 



{(p * n){z) := / (p{z - y)n{dy) < Ci(p{z), 
Jo 

for some < Ci < oo. Then for any S > 



z>0 



Sect{S^) < 1 + 



2Ci 


dn 




a 


dn 


oo 



(3.15) 



(3.16) 



Proof. We may assume that ||(in/(in||oo < oo. By virtue of Theorem 3.3, it is enough 
to show (3.10) with C = Cir for any r > such that n{dy) < rn{dy). Applying 
Schwarz's inequality, we can dominate BQ{f,gY by 

j A{dx)x j n{dy)f\x + yf ■ j A{dx)x j n{dy){g{x + y) - g{x)f 

< j dx(p{x) J n{dy)f'{x + yf ■ r j A{dx)x j n{dy){g{x + y) - g{x)f . 

Since by (3.15) 

j dx^ix) J n{dy)f'{x + yf ^ J dzf'{zf{ip * n){z) 

< C,j dzf'{z)Mz) = Ci y A{dz)zf'{zy, 

the desired inequality is derived. ■ 

A key ingredient to verify (3.15) is the following fact taken from Eq. (6) in the 
proof of Lemma 2.6 in [21]. (The function K there is identical with (p in the present 
paper.) 

Lemma 3.5 In addition to the assumptions in Theorem 3.3, assume that a > and 
0<&:=6 + c<oo. Then the Levy measure A in (3.3) has a strictly positive density 
dkjdy —: p>{y)/y with </? being dijjerentiable. Moreover, (f{0) :— limyi^o (f{y) — 1/a 
and 

{(p*n){y) = a(f'{y) + b(p{y), y > 0. (3.17) 

To grasp the validity of (3.17), it is worth noting that taking the Laplace transform 
of both side of (3.17) leads, at least at formal level, to the equation equivalent to the 
one derived by differentiating (3.3) with q = provided that ip{0) = 1/a. Note also 
that (p'{0) :— liuhyi^o (p' (y) — —b/a^ is deduced from (3.17). 

Proposition 3.6 Under the same assumptions and with the same notation as in 
Lemma 3.5, define V{y) — — \og(p{y) (, so that\\va.yy^V'{y) —: V'{0) exists). If 

sup{V'{0) - V'{y) : y > 0} < C2 

for some < C2 < 00, then for any 6 > 



Sect{£^) < 1 + 



\ 



2C, 



dn 
dn 



(3.18) 
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Proof. By letting 7/ ; in (3.17) a^'{0) + b^{0) = or 6 = -a^'{0)/^{0) = aV'{0). 
So again by (3.17) 

* n){y) = aip'iy) + aip{y)V'{Q) = 0(^(7/) (-l^'(t/) + ^'(0)) < aC^ifiy), 

and thus (3.15) with Ci = 06*2 holds true. Therefore (3.18) follows from (3.16). ■ 

In the reversible case n = 0, the function V in Proposition 3.6 is affine, so that we 
can take C2 = 0. Under some integrability condition on n, quantitative information 
of C2 can be obtained from Lemma 2.6 in [21] combined with Eq. (6) there. It will 
turn out that (3.18) is one of basic tools in the next section, where more specific cases 
are discussed. 

A naive guess based on the integration by parts formula (3.5) would be that 
Sect(£) = 00 whenever a = 0. We have not proved this, nor given any sufficient 
condition for the CBCI-process not to satisfy the strong sector condition. We just 
present a simple (but very special) example of such a CBCI-process. 

Example 3.1 Let 6, c > and n{dy) = cei{dy), where ei is the delta distribution 
at 1. Then, for each 5 > 0, the CBCI-process with quadruplet (0, 6, n, 5) is ergodic 

and does not satisfy the strong sector condition. Indeed, letting f{x) = sin27ra; and 
g{x) = cos27ra;, one can observe from (3.2) and (3.5) that S^{f,f) = = £^{g,g) 
and that 

^^{f^9) — "^T^clx ( cos2n(x + y)(cos2TT(x + y) — cos2TTx)dy)s 
Jo 

= 'Kc{x)s = 7rc(5$'(0) = TTcS/b, 

respectively. 

For later use, we close this section by giving a lower bound of Sect(£^) in a general 
setting as a refinement of the calculation (3.1). 

Proposition 3.7 Let S be as in Proposition 3.1 and f,g E D(L) be such that 
^(/, m{9: 9) > and S{f, g) := (£(/, g) - S{g, f))/2 > 0. Then 



Sect(£:)2 - 1 > 



^ (^(/,/)^(^,^)y^^(/,^)^(/,^),A(/,^)>0) 
00 ' {£{fJ)£{g,g)^£{f,g)£{f,g)A{L9)-^) (3.19) 

^^^'^^ {£{f,m9,9) = £{f,9m,9)), 



y s{f,9) 

where A(/, g) = £{/, f)£{g, g) - £{f, gf > 0. 
Proof. The proof will be based on 
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Put U(t) = S{f,f + tgY/£{f + tg,f + tg). By direct calculations it can be shown 
that ^ log U [t) vanishes only for 

t — If) 



S{f,mg,g)-S{f,g)S{f,g) 

if S{f,f)S{g,g) — S{f,g)S{f,g) ^ 0, whereas U'{t) never vanishes unless U{t) — 
if S{f, f)S{g, g) — S{f,g)£{f,g) = 0. In the former case, the supremum in (3.20) 
is achieved at t = to ^-nd (3.19) is obtained for the first two cases by calculating 
limt^to U {t). In the latter case, noting that £{f, g) = £{f, g) + £(/, g), we have 

U{t) y U{t) £{f,gy 



tenSifJ) l*Koo £(/,/) £(f,f)S{g,g) 



S{f,g)S{f,g) S{f,gy 

which proves (3.19) for the third case. Lastly, the inequahty A{f,g) > holds in 
general since £ is symmetric and nonnegative definite. ■ 



4 Generalized gamma convolutions as stationary 
distributions 

In this section, wc apply the results in the previous section to the sector constant 
estimate for a class of CBCI-processes whose stationary distributions are GGC's. 
(The general reference for GGC's is [2]. The interested reader is referred also to [22] 
or [6].) The situation, however, is a converse of that in Section 3 in the following 
sense. We shall be given a priori g > and the Levy measure A of some GGC, and 
then intend to choose a, b and n so that (3.3) holds. To be more specific, recall that 
a GGC is an infinitely divisible distribution on R_|_ with Levy measure of the form 

Amidy) := (| e-"^m(rf«)) ^, (4.1) 

where m (referred to as the Thorin measure) is a measure on (0, oo) with 
/ \ log u\m{du) + / m{du) < oo. 

J{0,l/2] .^{1/2,00) 

This condition is necessary and sufficient for the Laplace exponent 

*,,m(A) := gA + 1 (l - e-^^') Am{dy) =<1^ + j log (l + ^) ^{du) (4.2) 

to be finite for all A > 0. (Notice that found in the literature is the condition that 
/(o,i] \ logu\m{du) + J^]^^^^u~^m{du) < oo, which allows m{du) — \{Q^i){u)du/\\ogu\ 
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inappropriately. Here and in what follows, the notation 1e denotes the indicator 
function of a set E.) We call the distribution having the Laplace exponent (4.2) the 
GGC with pair {q, m). It must be remembered that every gamma distribution has a 
degenerate Thorin measure. 

In order to study the above mentioned problem, we do a heuristic calculation; 
differentiate (3.3) with Levy measure (4.1) to get 

r I 1 
q + / m{du) = , A > 0. (4.3) 

^ + '^ aX + b+ (1- e-^'')n(du) 

This equation motivates us to exploit the theory of Bernstein functions [18]. Defining 
Ai to be the totality of measures m on (0, oo) such that /(I + u)~^m{du) < oo, we 
recall that every complete Bernstein function g is represented uniquely in the form 

f A 

g{X)^qX + r + J j-^m{du), A>0 (4.4) 

where q,r > and m E Ai (cf. Remark 6.4 in [18]), and Proposition 7.1 in [18] asserts 
that a function g : (0, oo) — )■ R is a non-zero complete Bernstein function if and only 
if g*{\) := A/5'(A) is a complete Bernstein function. With the help of these two facts 
one can show the next lemma, a key to our construction of the desired CBCI-process. 
In what follows we adopt the convention that l/oo = and set = / u°'m{du) for 
q; e R. 

Lemma 4.1 Let q,r > and suppose that m E M. is non-zero. Then there exist 
uniquely a,b> and M e such that 

q + ^ + I —^m{du) = , A>0. (4.5) 



A J \ + u ^ ' ^ , I. , f ^ 



M(du) 



Moreover, 



(g > 0) h-[^ (r > 0) 

l/(r + mo) (g = 0), ~1l/(g + m_i) (r = 0) 



(4.6) 



and 



_ ' l/q-h (g > 0) 

Mo = < mi/(r + mo)^-6 (g = 0,mo<oo) (4.7) 
oo (gr = 0, mo = oo). 

Proof. The first half is immediate from the general facts on Bernstein functions 
previously mentioned. Indeed, defining the function g on (0, oo) by (4.4), we can find 
uniquely a, 6 > and M & M. such that 

^=,*(A)=aA + 6 + /^M(<i«), A>0. 
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This is nothing but (4.5). Most calculations needed to show (4.6) and (4.7) are simple. 
While letting A — )■ oo in (4.5) yields a — whenever g > 0, letting A | in (4.5) 
gives the value of b in (4.6). In the case q ~ 0, letting A — >■ oo in (4.5) multiplied by 
A shows that a = l/(r + mo)- With (4.6) in mind (4.7) can be proved in a similar 
manner because Mq = lim^^oo/ j^M^{du). For instance, in the case where q — 
and mo < oo, (4.5) and (4.6) together yield 



Letting A — > oo proves (4.7) in this case. The proof for the other cases are left to the 



Denote by S{m) the support of a measure m. (4.7) together with Schwarz's inequality 
implies that M = if and only if g = = r and m is degenerate. In such a case, 
it is understood that S{M) = 0. The next lemma gives bounds of inf »S(M) and 
sup5(M) in terms of q, r and m. 

Lemma 4.2 Let q,r > and suppose that m & M. is non-zero. Let a,b > and 
M & A4 be as in Lemma 4-i- Define s_ = •§-(?, r, m) and s+ = «+(?, m) by 




reader. 




and 




respectively. Then the following assertions hold, 
(i) If a — — b, then 



r 



inf cS(m) < s_ < inf 5(M) 



r + mo + q inf S{m) 



and 




(a) If a > and b — 0, then 



r 



inf 5(m) < s_ < inf 5(M) < sup5(M) < sup5(m). 



r + mo 
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(Hi) If a = and b > 0, then 

miS{m) < miS{M) < supc5(M) < 



< sup S{m) + moq ^ > 0) 



(iv) If a,b > 0, then S{M) is contained in the closed interval from inf»S(m) to 
sup<S(m). 

Before going to the proof, it is worth noting that the analytic extension of (4.5) can 
be regarded as a relation between Stieltjes transforms of m and M. Introducing the 
notation Gm{z) — J{z — u)~^m{du) for m e A4, we deduce from (4.5) 

Gmiz) - g + - = , ^ ^ , ^ e C \ [inf 5(m),sup5(m)]. (4.8) 

Z OjZ — — ZKjrM\Z) 

Proof of Lemma 4-2. We shall employ the following fact. (See e.g. Theorem A. 6 
in [10].) For any si < S2, a measure M is supported on [si,S2] if and only if Gm 
is holomorphic on C \ [si, S2], negative on the interval (—00, Si) and positive on the 
interval (s2, 00). 

To show (i) we assume, in addition to a = 6 = 0, that < inf 5(m) < sup S{m) < 
00. By (4.8) we have 

Gm{z) = ^ 



/ f m(du)\ ( f m(du)\' 

U-/^^ -r r-z{q+^^] 
\ J z — u J \ J u — Z J 



from which it is easily verified that Gm is holomorphic on C \ [s_,s+], positive for 
all z > Sj^ and negative for all z < S-. So, s_ < m.iS{M) < sup (S(M) < s+ for 
the abovementioned reason. To prove the estimates for s_ and s+, we may assume 
further that mo < 00. We then have g > because of (4.6). It is clear that a 
lower estimate for s_ is given by the smallest solution, say s'_ — s'_{q,r,m), to the 
quadratic equation in s 

s( + ^0 \ ^ ^ 
\ inf (S(m) — s J ' 

namely 

s_ > s'_ 

(r + mo + Q inf S{m)) — ^J{r -\- mo + ?inf S{m)Y — 4gr inf S{m) 



2q 



(4.9) 



— = TTTT — T i^^^ S(m) . 

r + mo + q ml S[m) 



An upper estimate for s+ is given by the largest solution, say s', — s',{q, r, m), to 



s\q + 



sup<S(m) — 
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Thus s+ is dominated by 

(r + mo + ? sup <S(m)) + \J{r -\- rf%o + q sup <S(m))^ — Aqr sup <S(m) 



< sup <S(m) + 



2g 

r + mo 



Also, it is obvious that s+(= inf 0) = oo for g = 0. 

Next we prove (ii), assuming that a > and 6 = 0. By the former q = and 
mo < oo. Therefore, again by (4.8) 



Gm{z) = 



^ m{du) 



z — u 



, m(du)\ , . 

r - z I — — — I [r + mo) 



This allows proceeding along the same lines as the proof of (i). The proof of (iii) 
follows very closely the proof of (ii). So, these proofs are omitted. It remains to prove 
(iv). To this end, assume that a, 6 > and observe from (4.6) and (4.5) that 

m_i / — ^ — midu) — mo / , — rmidu) 
J z — u J uiz — u) 

Gm{z) = 



mom_i / m{du) 



u 



in which both mo and m_i are positive and finite. This shows the analyticity of Gm 
on C\ [inf S (m) , sup S (m)] . We may assume additionally that m is nondegenerate, for 
otherwise the assertion is trivial since M = 0. Then the positivity of the numerator 
on the right side for any 2; G R \ [inf S{m), sup S{m)] follows from 

m_i / midu) — mo / —, rmidu) 

J z — u J u[z — u) 

= I I — V r~ — AT^{du)m{dv) 

J J \v[z — u) v{z — V) J 

— [ midu) [ —. rmidv) — ( f midu)] > 0. 

J z — u J v{z — v) \J z — u J 

These properties together prove (iv). ■ 

We now present the main result of this section, which concerns the construction of 
the CBCI-process having a given GGC as a stationary distribution. In other words, 
the problem (II) addressed in Introduction is solved. Simultaneously, the sector 
constant estimate will be obtained as a solution to the problem (III) . Note that every 
Thorin measure belongs to M.. 

Theorem 4.3 Let g > and suppose that m is a non-zero Thorin measure. Let 
a,b> and M & M. he as in Lemma 4-1 with r — 0. Then the GGC with pair {q, m) 



21 



is a unique stationary distribution of the CBCI-process with quadruplet {a,b,n,l), 
where n is a measure on (0, oo) defined to be 



n{dy) ^dy j u^e-''m{du). 



(4.10) 



If, in addition, q — 0,mi < oo and miS{m) > 0, then toq < oo, miS{M) > and 
the CBCI-process with quadruplet {a, b, n, S) satisfies 



Sect{£^) - 1 < ^(^^ - miS{m] 
for any S > 0. 



< 



mi 



miS{M) ^ \mo-miS{m) 



- 1 



(4.11) 



Proof Firstly, we claim / mm{y , y'^}n{dy) < oo. To see this, note that by (4.10) and 
Pubini's theorem 

j min{y, y^}n{dy) = j M{du)u^ j min{y, y^je-'^ydy 
— j M{du) j yv[)m.{l,y/u}e~^dy 

// ru roo \ 

M{du) i^J {y^/u)e-ydy + J ye-Myj . (4.12) 

It is elementary to verify that the integrand in the last expression is bounded above 
(and below) by a positive constant times 1/(1 + m). Therefore it follows from M e 
M. that / mm.{y,y'^}n{dy) is finite. Applying Pubini's theorem, we see that /(I — 
e~^y)n{dy) = A /(A + ?/)-^M(dii) for A > 0. Plugging this into (4.5) with r = yields 
(4.3). By integrating it 

+ J log ^1 + ~^ m{du) = J 



du 



au + b + j {I - e-''y)n{dy) 



A > 0. 



In view of (3.3) and (4.2), this proves the first half. 

For the proof of the last half, we assume additionally that 9 = 0, mi < oo and 
inf S{m) > 0. The last two conditions together imply that mo, ^-i < oo since m is a 
Thorin measure. Hence a, fe > by (4.6). Also, by (4.10) and (4.7) / yn{dy) = Mq = 
mi/mp — l/m_i, which is finite. We now apply Proposition 3.6 to the Levy measure 
(4.1) or equivalently to the function (fiy) := / e~'^^m{du). It is easily observed that 
V{y) = — log 99(1/) satisfies 



V\0) - V\y) = 



J um{du) j e '^^um{du) 



m{du) 



-uy 



m{du) 



< ^ — inf S(m) 
mo 



for all y > 0. On the other hand, it follows from (4.10) that 



dn 



dn 



< 



< 



1 



inf5(M) - inf«S(m)' 
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where the last inequahty is imphed by Lemma 4.2 (iv). Therefore, (4.11) is deduced 
from (3.18). ■ 

Notice that the upper estimates (4.11) are effective in the sense that the most right 
side vanishes only in the reversible case, namely when m is degenerate. The symmetry 
found in the statement of Lemma 4.1 allows one to show a converse of Theorem 
4.3. Roughly speaking, we will see below that every ergodic CBCI-process with n 
having a (non-zero) completely monotone density has some GGC as a (non-reversible) 
stationary distribution. 

Theorem 4.4 Let a,b > and suppose that a non-zero M ^ M. is given. Define n 
and $ by (4-iO) and the first equality in (3.3), respectively. Then J mm{y , y^}n{dy) 
is finite and the following assertions hold true, 
(i) $(1) < oo if and only if b > or 

M{dv)Y\ 
— ^ — - du < oo. 
u-\- V J 

(a) //"$(!) < OO; then the unique stationary distribution of the CBCI-process with 
quadruplet (a, b, n, 1) is a GGC with some pair {q, m) such that 





_ (a > 0) _ _ 
1/(6 + Mo) (a = 0), """-^ 



' l/a _ (a > 0)_ 
Mi/(6 + MqY {a = 0, Mo < oo) (4.13) 
oo (a = 0, Mo = oo). 



(Hi) If a,b > 0, Mo < oo and infiS(M) > 0, then for any 6 > the CBCI-process 
with quadruplet {a,b,n,6) is ergodic and (Sect{£^) — ij is dominated by 



(a inf S{M) - 6 - Mq) + 4aMo inf S{M) - (a inf S{M) - 6 - Mo) 



(4.14) 



ainf5(M) 

Proof. By virtue of (4.12), that / min{y , y'^}n{dy) < oo is ensured by M G Ai. 
Putting guiu) = f{u + v)^^M{dv), note that $(A) = f^i^au + b + ugM{u))^^du. For 
the proof of (i), it is sufficient to show that 6 = and fl{ugM{u))~^du = oo together 
imply $(1) = oo. If M_i < oo, this holds true since $(1) > jQiau+n'M-iY'^du = oo. 
If M_i = oo, there exists v e (0, 1) such that a < gniu) for any u e (0, v] and hence 
$(1) > Jq {2ugM{u))~^du = oo. These observations prove (i). 

To show (ii) we interchange in Lemma 4.1 the roles of (a, 6, M) and (g, r, m) to 

get 

a + \+ j -^M{du) = ^ , A>0 (4.15) 

qX-\-r-\- / midu) 

J X + u ^ ' 

for some g, r > and m € }A. Here, according to (4.6), q is given by (4.13), r = 
if 6 > 0, and r = l/(a M_i) if 6 = 0. But it follows from $(1) < oo and (i) that 
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M_i = oo whenever 6 = 0. As a result r = 0, and therefore the formula (4.13) for 
mo precisely corresponds to (4.7). By (4.15) and Pubini's theorem one can show that 
$(A) = ^q^rni^)- Since $(1) < oo, this proves not only that m is a Thorin measure 
but also the assertion (ii). 

It remains to show the sector constant estimate under the stronger assumptions 
that a, 6 > 0, Mq < oo and miS{M) > 0. This should be derived from (4.11) once 
nil/ mo and inf S{m) are estimated in terms of a, h and M. For this purpose, combine 
mo = l/o with Mq — rni/fno^ — b to get mi/mo = (-^o + b)/a. Also, Lemma 4.2 (i) 
and (4.9) with (g, r, M) being interchanged with (a, b, m) are applied to derive 

mlS{M) > s4a,b,M) > s'_{a,b,M) 

(a inf S{M) + b + Mq) - ^(ainf 5(M) + 6 + Mq)^ - 4a6inf 5(M) 
~ 2a ■ 

These calculations together yield 

^ — inf 5(m) 
mo 

{b + Mo) - a inf S{M) + ^ (a inf 5(M) +b + Mq) ^ - 4.ab inf S{M) 
~ 2a 

^(ainfcS(M) - 6-Mo)%4aMoinfcS(M) - (a inf cS(M) - 6 - Mo) 

~ 2a ■ 

Therefore, the bound (4.14) is deduced from the first inequahty in (4.11). ■ 

To check, consider the case discussed in Example 2.1 (ii). It is easily seen that 
the Laplace exponent (2.23) is that of the GGC with q — {) and Thorin measure 

m(du) = — — i . • (4.16) 

^ ' r(Q;)r(l-Q;) {u-kY ^ ' 

Accordingly, mo = oo and m_i = So, (4.6) gives a = and b = consistently. 
Although Lemma 4.1 itself does not give any explicit form of M, we can identify it 
with 

^(*' = f9^ 

thanks to (2.22). Therefore, in this case S{M) = S{m) = (/t, oo). In the next section 
we will be provided with some procedure to derive a formulae for M including (4.17) 
via (4.5) and with further examples as well. 

As for the lower bound of the sector constant of CBCI-processes discussed in 
Theorem 4.3, one can show the next result as an application of Proposition 3.7. 

Theorem 4.5 Suppose that m is a non-zero measure on (0, oo) with inf <S(m) > 
and mo < oo. Let a,b > and M & M. be as in Lemma 4-1 with q — — r and n 
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be given by (4-10). Then the bilinear form £ associated with the CBCI-process with 
quadruplet (a, b, n, 1) satisfies 

fm_im_3 - rrf_2) 

Sect(£) - 1 > - — . 

~ 2m_imi2^-3 + 4mi^ml2 + 12m^;^m_2m_4 — Qmi^^m^g — 

Proof By the assumption m is a Thorin measure and the GGC with pair (0,m) has 
an exponential integrabihty. Indeed, its Laplace exponent $o,m(A) can extend real 
analj^ically to A > — inf S{m). In particular, we have the finite moments given by 

{x') = (-1)^ S-.e-''^'-^'^ 

A=0 

for A; = 1, 2, ... . Since $5*1(0) = {-lf-^{k - l)!m_fc, it follows that 

{x) = m_i, (a;^) = m\ + m_2 and {x^) = 2m_3 + 3m_im_2 + nf^i- (4-18) 

We are going to apply (3.19) by taking f[x) = x and g{x) = x"^, for which A(/, g) > 
holds by virtue of (3.2) and Schwarz's inequality together. Explicitly, our task is to 
show for some common C > that 

E{f, gf = C(m_im_3 - m\f (4.19) 

and that A(/, g) = S (/, f)S {g, g) - S (/, gf equals 

C(2m_im^2^-3 + + 12m^im_2m_4 — dTfTtirfT^^ — m^2)- (4.20) 

Noting that / n{dy)y^ = k\Mi_k {k = 1,2,.. .) by (4.10), observe from (an extension 
of) (3.5) and (3.2) that 

^{Lg) = \{x j n(dy)(-y')) = -^{x) J n{dy)y' = -(x)M_2, (4.21) 
S(f, f) = (^) (« + ^ / n(dy)y'^ = {x) (a + M_i) , 

£{f,g) = 2{x') {a + ^J n{dy)y'^ + ^{x) J n{dy)y' 
= 2(x2) (a + M_i) +3(x)M_2 

and 

£{g,g) = A{x')[a + ^Jnidy)y'^+2{x') Jnidy)y'+^{x) Jn{dy)y' 
= A{x^) (a + M_i) + 12(a;2)M_2 + 12(a;)M_3. 
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The last three equaUties together yield 

+12{xf (a + M_i) M_3 - 9(x)^M_2^ (4.22) 
We shall calculate M_i, M_2 and M_3. The results are 

mom_i m_i 



and 



m 



-1 

whose proof are postponed for a while. Then (4.19) with C = l/m^^ follows imme- 
diately by plugging (4.18) and (4.23) into (4.21). In principle, (4.20) with the same 
C can be obtained similarly though the calculation is tedious. To carry this out, 
observe that a + M_i — rn-2/rffii by a = 1/mo and that (4.18) give 



{x){x ) — {x ) = m_^m^2 + 2m_im_3 — m_2- 

(4.20) with C = 1 /m^i will now be derived in a fairly straightforward way. 

It remains to prove (4.23) and (4.24). We exploit a variant of the identity used 
to find Mq in the proof of Lemma 4.1 (with q — — r): 



r 1 / T— midu 

/ -J—M(du) = J-A + u 



_ /■ A , , m_iA' 

mo / m[au) 

J X + u 



^ A > 0. (4.25) 



X + u 

Letting A | 0, we get (4.23) for M_i with the help of L'Hospital's rule. By differen- 
tiating (4.25) 



um{du) f m{dv) _^ f f 'm{,du) 



1 J {X + uy J V \J X + u 



(A + «)2 ^ ' A2 _ f f m{du) 

m_i 



X + u 



II 



(11^ — y\'^ 

jm{du)m{dv) 



{X + u)'^uv{X + v) ^^26) 



, , midu) 



X + u 



where the symmetry of m{du)m{dv) applies to show the last equality. Letting A 4- 
leads to (4.23) for M_2. Finally, differentiating the both sides of (4.26) multiplied by 
the square of J{X + u)~^m{du) and then letting A 4- 0, we see without difficulty that 

ir^ 9 ^-irr 2(m-im^A — m_2m_3) 

2M_3m^ + 2M_2m_im_2 = ^ — ^ -. 

m_i 

This combined with (4.23) proves (4.24). The proof of Theorem 4.5 is complete. ■ 
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5 Further discussions and related topics 



In this section, most calculations will be based on the equations (4.5) and (4.8) 
with r = 0, and therefore we take r = without explicit mention (except in the 
statements of Propositions 5.2 and 5.3 below). The first two subsections are devoted 
to further studies of the correspondence between (a, b, M) and (g, m) under some 
special circumstances. In the final subsection, the basic relation (4.5) will be discussed 
also in connection with certain topics in noncommutative probability theory. 



5.1 Discrete Thorin measures 

This subsection concerns the correspondence between GGC's and ergodic CBCI pro- 
cesses under the condition that m (or M) is discrete. As the simplest example of 
such GGC's (other than gamma distributions) we first discuss the case of convolu- 
tions of two gamma distributions in rather detail. Given a; G R, denote by the 
delta distribution at x. 

Example 5.1. (i) Let 71, 72, Ai, A2 be positive constants. Consider m :— ^yie^^ +^2^X2 
as the Thorin measure. According to (4.6) with q — 0, a and b are chosen as 

a — and = -— j — i — , (5.1) 

7i + 72 Ai 7i + A2 72 

respectively. Our ansatz here is that the measure M satisfying (4.5) with g = is of 
the form M — ce^, for some c, k > 0, which are to be determined. The measure n in 
(4.10) is then given by n{dy) — CK^e~'^'^dy and the equation (4.5) reads 

+T^ = A>0. (5.2) 



A + A, A + A2 ^^^^^ cA 



It is not difficult to see that the above requirement is fulfilled by setting 

a A271 + Ai72 , . 

K := -A1A2 = ■ (5.3) 

b 71 + 72 

and 

M ^\ ^ A 7i72(Ai - A2)^ , . 

c :— a Ai + Xo — K) — b— ttt-t ; -. 5.4 

^ ^ (7i + 72)'(A27i + Ai72)^ ^ ^ 

which vanishes for degenerate m with Ai = A2 in accordance with the comments in 
the paragraph preceding to Lemma 4.2. Lastly, assuming Ai < A2 and letting 5 > 0, 
we have by the first bound in (4.11) combined with (5.3) 



Sect{S^) - 1 < 



71A1+72A2 , \ 2 /2(A2-Ai)72 



^ V 7i + 72 J 1^ V A271 + Ai72 

for the CBCI- process with quadruplet (a, b, n, 5). 
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(ii) One can reverse the above procedure to construct explicitly the Levy density of 
the stationary distribution of the CBCI-process with n being of the form n{dy) = 
CK'^e~'^'^dy. Indeed, given a,b,c,K > 0, we find Ai,A2, 71,72 satisfying (5.2) in the 
following manner. In view of the first equalities of (5.3) and (5.4), the required Aj's 
must solve the equation p{X) :~ aX^ — {an + 6 + c)A + 6/« = 0. This leads to 

iaK + h + c) - VD , ^ {aK + h + c) + ^/D 

Ai = and A2 = , 

la la 

where D — {an -\-h-\- cf — AabK = {an — b — + AacK > 0. Note that Ai < k < A2 
since p{k) < 0. Moreover, 7j's are determined by 

K — Xi , A2 - K 

7i = -7T TT 72 



a(A2-Ai) a(A2-Ai)' 

for which (5.1) are easily checked to hold. Consequently, the Levy density (4.1) equals 
(716""^^^ + 726"'*'^^) /^/. In addition, one can derive from (5.5) the sector constant 
estimate described in terms of a, b, c and k by noting that 

2(A2 - Ai)72 _ 2(A2 - k) _ y/D -{an-b-c) 
A271 + Ai72 K an 

which clearly corresponds to (4.14). 

Example 5.2. Given g > and a degenerate Thorin measure m, we have a similar 
situation to Example 5.1 except the sector constant estimate, which is not available 
since a = by (4.6). Indeed, for m = 7ieAi with 71, Ai > being fixed arbitrarily, 
6=1/ ((7 + 71/Ai) by (4.6) and it is easily verified that the measure M satisfying (4.5) 
is given by 

M = i ■ T— J ex-,+^i/q. 

Conversely, if we are given 6 > and M — ce^ with c, k > 0, then the measure m 
determined by (4.5) with q — l/{b-\-c) and a = is shown to be 

CK 

Alternatively, this can be derived directly from (3.4) by noting that n*^{dy) — 
{cK)^y^-^e-''ydy/{N - 1)! for = 1, 2, ... . 

Apart from explicit expressions, the above examples are generalized as follows. 

Proposition 5.1 Let I > 1 be a fixed integer, 
(i) For every discrete measure 

= 7ieAi H with 7i > anc? < Ai < • • • < A;+i, (5.6) 

there exists a unique measure M of the form 

M — cie^j + • • • + cie^i with q > and < ki < ■ ■ ■ < ki (5.7) 
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satisfying (4-5) with a = l/mo, b = l/m_i and q = = r. Furthermore, the 
sequences {Aj} and {/tj} interlace: 

Xi < Ki < X2 < ■ ■ ■ < < Xi < Ki < Xi+i. (5.8) 

Conversely, for given a,b > and a discrete measure (5.7), there exists a unique 
measure m of the form (5.6) satisfying (4.5) with q — — r. 
(a) Let q > be given. For every discrete measure 

^ — Ti^Ai + • • • + Ti^Ai with 7j > and < Ai < • • • < A;, (5.9) 

there exists a unique measure M of the form (5.7) satisfying (4-5) with r — — a 
and b = l/{q + rn-i). Furthermore, it holds that 

Ai < < A2 < • • • < Ki-i < Xi < K,i < Xi + moq~^. (5.10) 

Conversely, for given b > and a discrete measure (5.7), there exists a unique 
measure m of the form (5.9) satisfying (4-5) with q — l/{b + Mq) and r — — a. 

Proof, (i) Observe from (4.5) with q = that g{X) := J {X + u)^^ M (du) is a rational 
function of the form 

where a and b are given by (4.6) and 

Z+l I 

P(A) = + A,), g(A) = E 7^ + A,). 

i=l 1=1 jy^i 

We see also that Po{X) := (-P(A) — (aA + 6)Q(A))/A is in fact a polynomial with degree 
less than or equal to / — 1, and that Q has a zero in the interval (— Aj+i, —Xi) for 
each i e {1, . . . , /} because Q{—Xi)/Q{—Xi^i) < as observed from (5.6). Therefore, 
g{X) = aPo(A) /((A + Ki) • • • (A + ki)) for some Ki, . . . ,ki satisfying (5.8). It remains 
to find ci, . . . , Q > such that 

I 

aPo(A) = E Q n(A + Kj), A > 0. 

Since ki, . . . , k/ arc mutually different, a necessary and sufficient condition for the 
above identity to hold is that 

aPo{-Ki) ^ Ci'[[{-Ki + Kj), ie{l,...,l}, (5.11) 

which uniquely determine ci, . . . , q. Noting that Po{—Ki) — P{—Ki)/{—Ki) because 
of Q{—Ki) — 0, we can verify the positivity of q's by making use of (5.11) and (5.8). 
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It is almost a routine matter to show the converse assertion, whose proof we shall 
sketch. Let M be given by (5.7). Defining this time 

; I 
P{X)^l[{X + Ki) and Q{X) = J2c,]l{X + Kj), 

i=l i=l j^i 

we only have to show that the rational function 

can be rewritten into Z!-=} 7?/ (A+Aj) for some 7i > and Aj > satisfying (5.8). Such 
A2 . . . , A; are found as zeros of Qq since (5o(— '«i)/Qo(— '^i+i) < {i — 1, ... ,1 — 1). 
Therefore 

go(A) = (A + A2)---(A + A0gi(A) (5.12) 

for some quadratic polynomial Qi(A) = aA^ + giA + q2 with gi G R and q2 > 0. 
Moreover, with the help of (5.12) and (5.8), we can show that Qi{—Ki) < for each 
i e {1, . . . , These observations imply that Qi{X) = a(A + Ai)(A + A;) for some 
Ai e {0,K,i) and A/+i e («;;,oo). The rest of the proof (i.e. finding 7j's) is the same 
as before and the details are left to the reader. 

The assertion (ii) can be shown in an analogous way to the assertion (i). So we 
omit the proof of (ii) except the last inequality in (5.10), which follows immediately 
from Lemma 4.2 (iii). ■ 

5.2 Absolutely continuous Thorin measures 

As for continuous m, combining (4.8) with the following inversion formula of Stieltjes- 
Perron (cf. [24], p. 340, Corollary 7a) may provide explicit information on M: 

^^ M({s})+M({t}) 1,. fK ^ , 

M{(s,t))-\ \ ii-i^ = hm / lmGM{x + iy)dx, -00 <s<t< 00, 

2 TT 3/4-0 Js 

where i = \/— 1. In particular, M has a density function given by 

limlmGM(a:^ + ly), x e [s,t\, 



dx 71 yio 

provided that the right side converges boundedly and pointwise on [s,t]. We try to 
evaluate the right side in our setting. Let m e M. he non-zero and set q — for 
simplicity. Rewrite (4.8) in the form Gm{z) — —l/{zGm{z)) + a — b/z, where a, 6 > 
and M ^ Ai correspond to the pair (0,m) in the sense of Theorem 4.3. In the case 
where m is absolutely continuous, it would be possible to establish a density formula 
for M under suitable conditions, which do not seem, however, to be prescribed in a 
neat fashion. In typical cases, there exists a signed measure m' on [0, 00) such that 
Gm{z) — /[o,oo) log(2; — u)m'{du) and hence 

\im.Gjn{x -\-iy) = / \og\x — u\m' idu) + mm! iix 00)) 
j/4-0 J [0,00) 
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for each x > with m'{{x}) = 0. Therefore, the density formula would take the form 



dM 1 , 1 

lim Im 



dx nx yio Gm{x + iy) 

1 m'{{x,oo)) 

^ ~'Tf V 2 ^ ^ ^ 

/ log \x — ulm'idu) ] + (nm'Ux^oo))) 

\J[0,oc) J 

for X > 0. Similarly, (4.8) in principle makes it possible to derive the information 
of m corresponding to given a,b > and M E Ai. We shall give some concrete 
examples of such a kind, in which we continue to take q = 0. 

Example 5.3 (i) For < o; < 1 and k > 0, let m be as in (4.16). Then Gm{z) = —{k — 
z)~". Here the power function 2;" of 2; G C \ (— oo, 0] is defined to be |z|"(cosarg^ + 
isinarg^;) with arg^; chosen so that | arg2;| < vr. Since we know from Lemma 4.2 (iii) 
that S{M) C [k, 00), fix an a; > K arbitrarily. It is easy to see that 

limlm — — — = (x — k)" sin(Q;(— tt)) — '^"^ '^^ 



nx yio Gm{x + ly) nx x r(Q;)r(l — a) 

It is not difficult to verify that the above convergence is uniform in x on every compact 
interval contained in (k, 00). Thus (4.17) has been recovered. 

(ii) Let < Ai < A2 be arbitrary. Define m{du) = l(^\-^^^x^){u)du, for which a = 
l/(Ai-A2), b= l/(logA2-logAi). Moreover, Gjn{z) = log(2; - Ai) - log(2; - A2) and 
hence m' = — ex^. By (5.13) 

^ _ 1 . l(Ai,A2)(^) .5 

dx X / x-XA\ ^ 



While (4.7) tells that 



M([A„A2])= 



2(A2 - Ai) logA2-logAi' 

it is not clear how to verify this directly from (5.14) unless Ai = 0. If Ai > 0, we have 
also the sector constant estimate (4.11) which reads Scct(£'^) — 1 < \J (A2 — Ai) / Ai 
for each 5 > 0. For the special choice Ai = and A2 = 1, the density function of the 
GGC with pair (0,m) can be found in [6] (Eq.(259)). 

(iii) Let a, 6 > be given arbitrarily. Consider M{du) = l(^o^i){u)du for simplicity. 
The density of the absolutely continuous part rric of the Thorin measure m corre- 
sponding to (a, b. M) in the sense of Theorem 4.4 (ii) can be calculated from (4.8) with 
q = = r, namely Gm{z) = 1/ {az — b — zGm{z)), where Gm{z) = logz — log{z — 1). 
Indeed, defining H{z) — a — (b/z) — Gm{z) for -2 7^ 0, one can show that, for each 
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X > with X 1, H{x + iy) — )■ a — {b/x) — log|a;/(x — 1)| + i7rl(o,i)(a;) as y I 0. 
Therefore, 



drric 1 

— — = hm Im 

ax 



1 



1 



1(0,1) (a^) 



TTX J/4-0 H{x + iy) x / h x 

[a log 



X 



\ — X , 



+ 7r2 



The point which requires extra care is the unique pole, say Xq, of G^, which is located 
on the interval (1, oo). It is characterized as a unique solution to 



ax — h-\- X log ^1 ^ — 0, 



(5.15) 



from which h/a < xq < l + (6+l)/a can be deduced with the help of elementary 
inequalities u/{l + u) < log(l + u) < m for m > —1. The point mass of m at xq is 
given as the residue of Gm at z — xq' 



"^({^o}) = Urn Gm{z){z - Xq) 



z—>-xo 



(a + log f 1 - — ) + - 

V \ XnJ Xn 



d_ 

dz 

-1 



[az — h — zGm{z)) 



b 1 

+ 



Z=XO'' 

-1 



Xq/ Xq. 

where (5.15) has been used to get the last expression. 



Xo Xo-1, 



5.3 Connections with noncommutative probability theory 

In the previous subsection, the calculus of Stieltjes transforms played quite an im- 
portant role. So, it might be no surprise that observations we had made so far have 
some connections with noncommutative probability theory, another context in which 
the reciprocal of the Stieltjes transform serves as one of essential tools. We will be 
particularly concerned with the Boolean convolution and the free Poisson distribu- 
tions (known also as the Marchenko- Pasture laws). To give the definition of the 
former, we follow [19] and introduce the 'Boolean cumulant' Km{z) := z — 1/Gm{z) 
for a probability measure m on R. For our purpose the domain of such an operation 
shall be restricted to J\4i, the totality of probability measures on (0, oo). For any 
mi, 777,2 £ A^i, the Boolean convolution mi l±l 7772 of 7771 and 7772 is then defined to be 
an element of A4i such that 

{z), z eC \ R+. (5.16) 

Let t > be arbitrary. Following [3], one can define also the tih Boolean convolution 
power 777^* of 771 e A^i by 

Km^t{z)^t-Kmiz), zeC\R+. (5.17) 

It is a good exercise to verify from Lemma 4.1 with a — 1 that the requirements (5.16) 
and (5.17) determine uniquely such measures 7771 I±l7?72 and 777'*'*, respectively. (In fact. 
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for the verification, one needs to observe in Lemma 4.1 additionally that M_i = oo 
whenever 6 = 0. But this can be seen from (4.25).) The following proposition, which 
may have a number of variants, can be regarded essentially as a reformulation of this 
fact in the language of GGC's and the corresponding CBCI processes. 

Proposition 5.2 Suppose that probability measures mi andm2 on (0, oo) are Thorin 
measures. For each i G {1,2}, let the quadruplet (l,6j,nj,l) is the one determined 
from the pair (0,mi) by Theorem 4-3. Then the following assertions hold, 
(i) mi l±l m2 is a Thorin measure and the GGC with pair (0,mi l±) m2) is a unique 
stationary distribution of the CBCI-process with quadruplet (1, bi + 62, ni + n2, 1). 
(a) For each t > m^* is a Thorin measure and the GGC with pair (0,mf*) is a 
unique stationary distribution of the GBGI-process with quadruplet (1, ini, 1). 

Proof, (i) According to (4.10), are of the form ni{dy) = dy f u'^e~^^Mi{dy) for 
some Mi e M.. Observe from (4.5) with a — 1 that 

ff^^(du)Y' a = 6. + a/^>0, A>0. (5.18) 



\J X + u J J X + u 

Combining the above inequality with (5.16), we get 

{mi\£m2){du)Y^ _ / r mi(^y^ _ ^ 



X + u J \J X + u ^ 

or /(A + u)~^{mi l±l m2){du) < /(A + u)~^mi{du) for A > 0. By integrating with 
respect to the Lebesgue measure dX over [0, 1] and then applying Fubini's theorem 
J log{l + u~^){mi[t)m2){du) < J log{l + u''^)mi{du) < 00. Hence mi l±)m2 is a Thorin 
measure. The last half of the assertion follows from the equahties for i = 1 and i — 2 
in (5.18). Indeed, summing up them leads to 

(mi Wm2)(rfM)\"^ , > N , X f{Mi + M2){du) 



A = (61 + 62) + A y 



X + u J J X + u 



or 

(mi l±l m2){du) 



I 



X + u ^ , , u ^ , ^ f {Mi + M2){du) 



A + (61 + 62) + A y 



x + u 

for all A > 0. This is sufficient for the proof of (i). 

The proof of (ii) proceeds along the same lines as that of (i) on noting that by 
(5.17) and (5.18) 



/ 



mf(du} 1 max{f-', 1} 



J X + u J X + u 

for any A > 0. The details are omitted. 
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The final topic is related to the free probability theory. In that theory, the coun- 
terpart (in an appropriate sense) of the Poisson distribution exists and is called the 
free Poisson distribution. It is, by definition, of the form 

p (MA ._ / (1 - I^Mdu) + (3paAu)du (0 < /3 < 1) 

^a,0m - I p^^^(u)du (P > 1) 

for some a > and /? > 0, where 

1 



Notice that Pa^/s is a Thorin measure if and only if /3 > 1. The formula for the 
Stieltjes transform of this distribution is 



z + a(l-p)- J (z + a(l - p)y - Aaz 
Gp Jz) = ^ — ^- ^ — . (5.19) 

"'^^ ' 2az ^ ' 

(See e.g. p. 205 in [13] .) We now remark that a class of the free Poisson distributions 

plays a special role in describing the fixed points of the correspondence between m 
and M defined through Lemma 4.1 although we don't have any interpretation in the 
context of CBI-processes. Here is an explicit statement. 

Proposition 5.3 Let g > and suppose that m E M. is non-zero. Assume that 
a,b > and M E A4 satisfy (4-5) with r = 0. Then m coincides with M if and only 
if q = or mo < and m is given by 



m{du) 



' 1 — hq . , . , 

Pa,0[du) [mo < oo) 



« + g ' (5.20) 



— y'w - (6/2)21[((,/2)2,oo)(m)c^m {q = 0, mo = oo). 



where a — {1 — hq)/{a + qf and ^ = (1 + ah) /{I — hq). 

Proof. In view of (4.8), it is obvious that m = M if and only if Gm{z) — q 
l/{az — b — zGm{z)). By solving it we can deduce 



{a + q)z-b- J {{a + q)z - b)^ - 4(1 - bq)z 
Gm{z) = . 

In the case where mg < oo, we have a+q > by (4.6), and the proof concludes by com- 
paring this with (5.19). If mo = oo (and hence Mq = oo), then (4.7) and (4.6) imply 
q — Q and a = 0, respectively. Consequently, we have Gm{z) = {—b—\/b^ — Az)/{2z), 
from which the second expression in (5.20) can be derived by inversion. ■ 

Denoting by pq,a,h{u) and Ph{u) the densities of the measures on the right side of 
(5.20) in the first and the second cases, respectively, we note that 

limpo,a,6(M) = pb{u) = limp5,o,6(M) 

a\.\} qXU 
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for each m > and 6 > 0. Remark also that po{u)du coincides with m{du) in (4.16) 
with a = 1/2 and k = 0, for which clearly m = M holds. 

Acknowledgment. The author is grateful to Professors M. Yor and H. Masuda for 
bringing to his attention [18] and [8], respectively. 
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